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SUMMARY 



A fixed based simulator study of the approach and 
tangential landing of a class of manned re-entry gliders 
was conducted to determine the piloting techniques and ve- 
hicle parameters which best define the approach and land- 
ing problem. Tests were conducted on. 5 configurations 
typical of proposed vehicles with maximum lift to drag 
ratios of 3 to 5» The problem was Investigated from final 
approach to touchdown, In still and turbulent air, simulat- 
ing longitudinal equations of motion only. 

A three phase landing consisting of a steady glide 
final approach, a nearly constant load factor flare maneuver, 
and a 15 to 20 second shallow angle, final glide was the 
preferred technique. The vehicle parameters found to be 
most Important were maximum lift to drag ratio and speed 
for minimum drag. A landing analysis based on fitted para- 
bolic drng polorp ami a ©onolihnfe load faobor flare wa« 
found to be valid for the vehicles tested. 



INTRODUCTION 



The manned space vehicle capable of atmospheric re- 
entry and aerodynamic maneuvering to a horizontal landing 
is presently conceived as a fixed geometry wing-like lift- 
ing body of low aspect ratio and high profile drag. The 
vehicle may have relatively low wing loading and no pro- 
pulsion in the landing phase. 

The approach and landing of such vehicles is character- 
ized by high approach speeds, steep approach angles, rapid 
deceleration in near horizontal flight, and the probability 
of high speed and sink rate at touchdown. These character- 
istics pose a piloting task v/hich demands concentration 
and accuracy, while leaving scant margin for error. This 
is particularly true if landing facilities and navigation- 
guidance aids similar to those which accommodate modern 
high performance fighter aircraft are to be used. 

Approach and landing techniques which minimize the 
pilot judgment required, avoiding the usual flare or pull- 
up just prior to touchdown and providing reasonable capa- 
bility of pilot escape with today's ejection equipment, 
are proposed in REFS. 1, 2, 3* siid 5* technique has 

been used successfully with the X-15 research airplane 
(REF. 4). 
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The objective of this research was to study the tan- 
gential landing of a class of re-entry gliders using a 
fixed-base simulator with emphasis on determining the 
pilot and vehicle parameters which best define the approach 
and landing problem. 

Simulated landings of vehicles with maximum lift to 
drag ratios of 3 to 5 were conducted In still and turbulent 
air. The landing performance was analyzed, and an analytical 
procedure applicable to generalized vehicles developed. 

The principle pilot and author was a designated Naval 
Aviator with eight years of military flight experience and 
over 2400 hours In propeller type airplanes. 
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EQUIPMENT AND PROCEDURE 



1. Simulator. 

The simulator consisted basically of a fixed base seat, 
Instrument panel, CRT runway display, and side controller 
(FIG. 1). The computing equipment included an Electronic 
Associates Pace analog computer, a Llbratrol ^00 digital 
computer, with an analog-digital voltage converter. Data 
recording equipment Included a six channel Vlslcorder 906C 
strip recorder and the analog X-Y plotter. A detailed 
description of the simulator and associated equipment Is 
contained In Appendix A. 

The pilot was provided instrument display of the ve- 
hicle’s altitude, airspeed, pitch attitude, flight path 
angle, and glide path error. The runway simulation con- 
sisted of a porthole view of the runway boundaries shown 
In standard perspective. The simulated dimensions of the 
runway were 10,000 feet long by 200 feet wide. Airplane 
pitch altitude cues were represented by vertical motion 
and position of the runway display. Visual range and 
height cues could be Interpreted from the perspective 
view of the landing area. Only motion In the plane of 
symmetry was simulated and no lateral-directional task 
was assigned the pilot. 
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The side controller was given essentially constant 
force versus displacement characteristics with no centering 
forces. Breakout forces and dead zone were negligible. 

No control system dynamics were simulated. 

Atmospheric turbulence was simulated by introducing a 
rectangular or uniform probability density of angle of 
attack disturbances into the airplane equations. The 

disturbances could be converted to an equivalent vertical 



gust velocity by 






A vertical gust of 6.5 Tps RMS at the speed for minimum 
drag v/as used. 



2. Vehicle Dynamics. 

The basic longitudinal equations of motion for a rigid 
body with the standard axis system fixed in the aircraft 



are 


(REF. 6): 






(1) 






(Lift) 


(2) 


c/y 


(Drag) 


(3) 


/•^y — 1^3 0 




(Moment ) 




The kinematic equations are: 




(4) 


c 




(5) 












( 6 ) 



5 



where the terms sin ir' and cos ir have been linearized 
using a perturbation about , a nominal flight path angle. 



The aerodynamic forces 

~ f 

J)^ f (o<^ S') 



were expanded In a Taylor Series including second order 
terras where necessary to represent the re-entry vehicle 
aerodynamics. 




(8) 









o< y 



5i 






(9) 







ex' ^ 




:>cK^^ 



Aerodynamic forces were assumed Independent of Mach 
number for this experiment since the maximum Mach number 
simulated was approximately . 50 . Experiments with higher 
Mach numbers should consider compressibility effects. 

Equations 1 through 6 with 7» 8, and 9 were mechanized 
to simulate the vehicle dynamics In the experiment. A 
detailed description of the programming and the computa- 
tion of the runway display are given In Appendix A. 
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3. Vehicle Characteristics. 

In order to evaluate the possible effects of re-entry 
vehicle non-llnearltles , It was decided to fit as accurately 
as possible the wind tunnel data of typical vehicles, B and 
C (FIGS. 2 and 3). Tests were also conducted on a slightly 
simplified vehicle designated In which the second order 
aerodynamic coefficients were eliminated. Two additional 
hypothetical vehicles (A and D) of the same general class 
were tested. In each test the vehicle was programmed to 
simulate the desired characteristics In both a range of 
trimmed and unt rimmed flight conditions most appropriate 
to the Individual glider. 

The programmed characteristics of the five vehicles 
tested, A, B, B^, C, and D, are shown In TABLE I and FIG. 4. 
The characteristics for all vehicles tested are In the 
clean configuration. The lowering of landing gear and 
flaps was not simulated In this Investigation. 

4. Simulator Flight Task. 

The flight task was designed to simulate a straight 
In approach and landing on a 10,000 foot runway using 
essentially standard aircraft Instruments and visual dis- 
play. In preliminary experiments the problem was started 
with the vehicle at 20,000 feet altitude. However, the 
approach phase with the aircraft Initially positioned on 
a nominal approach path was found to be trivial and the 
test runs were shortened to save computer time. 
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Each test run was commenced In an equilibrium glide 
condition, Initially positioned on a nominal glide path 
to the runway at an altitude of 5*000 feet. The Initial 
velocity and flight path angle were varied over the range 
of practical flight conditions. 

In each approach configuration, the altitude at which 
the flare or pull-up was Initiated was varied until the 
landing performance for that configuration was considered 
optimum by the pilot. A series of at least 5 landings 
for each approach condition was then made In which the 
pilot attempted to execute the best possible approach, 
flare, and landing maneuver. 

The glider trajectory on all test runs was recorded 
by plotting V versus y on the X-Y plotter. Five of the 
following variables were normally recorded as a function 
of time: 







Cc/) 


A 




« 

r 


M 






X 


(f) 


0 



Pilot comments were made and recorded on each test run. 

No minimum speed for touchdown was specified. The 
maximum speed for touchdown was specified by an assumed 
maximum vertical speed of 10 fps. Landing performance 
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was Judged on controllability and accuracy for the flare 
maneuver, touchdown accuracy, ability to maintain vertical 
speed at touchdown within assumed limits, and overall 
consistency of the landing maneuver. Additional comments 
were made on the pilot's opinions of the short period ve- 
hicle dynamic characteristics where applicable. 

Simulated turbulence with an RMS equivalent to 6.5 fps 
for the vehicle average velocity during final glide was 
used in tests of vehicles A, B, and D. Pilot comments on 
the effects of turbulence were also made. 
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RESULTS AND DISCUSSION 

I. Description of the Landing Maneuver. 

An evaluation of the preliminary experiments validated 
the description .(REFS, 1 through 5) of the landing maneuver 
as divided into five parts or phases. 

The first part is a pilot guidance phase probably 
beginning after re-entry during which the pilot maneuvers 
to a high key position near or over the landing field aided 
by a radar-voice link and/or inertial navigation systems. 

The second part is usually described as a visual 
straight in or circling pattern in which the glider is 
maneuvered to Intercept the final approach leg at a pre- 
scribed altitude and airspeed. During this maneuver the 
pilot might use TACAN equipment as a primary navaid. 

The third and subsequent parts are the subject of 
this report and will hereafter be referred to as PHASES I, 

II, and III (FIG. 5). 

The three phases are described as follows: 

PHASE I - FINAL APPROACH. The final approach is a 
steady state glide at nearly constant airspeed and flight 
path angle aimed at a fixed, predetermined geographic point 
short of the runway. 

PHASE II - FLARE, The flare is essentially a constant 
load factor pull-up initiated at a predetermined altitude 
during which the sink rate is reduced to near the maximum 
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allowable for touchdown. The maneuver Is completed at a 
typically low altitude and small flight path angle. 

PHASE III - PINAL GLIDE, The aircraft Is flown In 
essentially cne "g”, constant ^ flight as excess speed 
and altitude are reduced for touchdown. During this phase 
the pilot has rather precise control over his altitude and 
flight path angle but essentially no control over decelera- 
tion. Landing gear and flaps are normally lowered In this 
phase. 

Another phase which Is not Included In this work Is 
the aircraft roll-out or slide-out after runway contact. 

2. Quality of Simulation. 

The preliminary evaluation of the simulator Indicated 
that the combination of Instrumentation and visual display 
provided sufficient realism to enable a pilot to perform 
and evaluate the final approach and landing maneuver. The 
absence of depth perception In the runway display was an 
early difficulty solved .by locating the flight path angle 
Indicator adjacent to the CRT display where It could be 
Included as a peripheral supplement to the limited physical 
depth cues. Accurate landings could then be performed with 
reference to the scope for attitude cues, and with airspeed, 
altimeter, and flight path Indicator providing precise con- 
trol of sink rate. 
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During the flare, the flight path indicator and altimeter 
provided the primary cues for controlling the flight path 
curvature in such a way as to arrive at the desired h and 
for final glide. The CRT display and attitude indicator 
provided the necessary attitude cues. 

3. Selection of Touchdovm Speed. 

Since the final approach and landing problem can be 
divided into 3 rather distinct phases, the analysis may 
be divided in the same manner by simply matching boundary 
conditions. It is also convenient to start the problem 
in reverse since it was found that the touchdown conditions 
can be determined independent of the approach, flare, and 
final glide. 

In general the selection of a minimum ground contact 
speed is not dependent on a conventional "stall speed". 

The re-entry vehicle, because it is most probably a lift- 
ing body, has no well defined Cj^ maximum. Instead the 
minimum touchdown speed might be determined by any one of 
the following: 

a. Wing buffet 

b. Directional stability limits 

c. Ground clearance limits 

d. Maximum landing gear extension 

e. Minimum trim speed 
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For this experiment It was assumed that there were 
no limitations on the minimum speed except the ability of 
the pilot to precisely control the touchdown. Elevator 
control power was provided In excess of that required to 
trim at the touchdown speed. Pitch damping was Incorpo- 
rated In all vehicles to Insure satisfactory short period 
dynamic characteristics at minimum speed. Since the long 
period vehicle dynamics are of paramount Interest In this 
study, difficulties with short period motions were thus 
virtually eliminated, (See Section 10 on Short Period 
Characteristics). Under these conditions the vehicle 
could be flown to a touchdown at any speed within the 
acceptable limits of vertical speed at touchdown assumed 
(10 fps). 

It was found In the experiment that If flight prog- 
ressed to a speed below a value characteristic of each 
vehicle, the loss of kinetic energy was so rapid that an 
excessive sink rate usually resulted. This Indicated that 
the pilot was unable to Increase the angle of attack rapidly 
enough to maintain one ”g” flight (L=W). This characteris- 
tic speed was found to be slightly below the speed for 
minimum drag (Vg). 

Assuming a characteristic parabolic drag polar where 

4 “ <4 



( 10 ) 



the for maximum lift to drag ratio Is 



( 11 ) 




Then 

( 12 ) 







and 

(13) 



ty ma'ii. 



/ 



Letting 

(14) 

(15) 






K 






=- 



A 



Substituting In 

(16) />= ^ 



yields 

(17) 








Then with the drag equation (2) 

(18) / 

jL 
Ji 






U' 



14 



In this equation the variable u represents a linear con- 
tribution of the profile drag to the deceleration while 
l/u is a hyperbolic contribution of the induced drag. 
The constant terms are a combination of gravitational 
acceleration and vehicle characteristics. 

For flight speeds less than Vg the induced drag pre- 
dominates and deceleration rapidly approaches infinity 
as showTi in FIG. 6. 

A simple statistical analysis of the successful 
simulator landings are shown in FIG. ?• Based on this 
empirical data a minimum touchdown speed for a reasonable 
probability of success should be approximately .88 Vg in 
the absence of other requirements for a higher landing 
speed as previously discussed. 



4. The Maximum Touchdown Speed. 

The maximum allov/able touchdown speed is generally 
a structural limitation imposed by the vehicle. The ver- 
tical speed at ground contact depends on both V and ^ 



since 



cJk 

dt ~ ^ 



FIG. 8 shoves h versus ^ where the limiting sink rate is 
assumed to be 10 fps. In practice the airplane will be 
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landed with some small but finite ^ found to be minus 2 
to 2.5 degrees in this experiment. PIG. 8 shows a dis- 
advantage of vehicles with high Vg, hence high landing 
speeds, requiring relatively smaller ^ at touchdown. 

It was found with all the vehicles flown in this 
experiment that a maximum sink rate of 10 fps is a prac- 
tical limitation and imposes no piloting problems unless 
landings are attempted at speeds below .88 Vg. These 
landings usually occurred unintentionally if the pilot 
leveled the vehicle high and subsequently ran out of kinetic 
energy before he could lower the vehicle to the ground. 

5 . Final Glide. 

The final glide (PHASE III) was found best character- 
ized by the time T between the completion of PHASE II and 
touchdown. This "tlme-to-go” to touchdown is the interval 
needed by the pilot to correct the errors in Judgment and/or 
execution which occur during flare allowing him to pre- 
cisely lower the rapidly decelerating glider to the runway 
with an acceptable sink rate. 

Experience with the X-I 5 airplane (REF. 4) and flight 
simulator experiments with fighter aircraft in high drag 
configurations (REFS. 1, 2, and 3)» indicates that a min- 
imum of 15 to 20 seconds should be allowed in the final 
glide. Reference 3 states that 30 seconds is usually the 
maximum required for any condition. 
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In the simulator flights, 10 seconds or less was found 
generally unsatisfactory resulting In unacceptable "scatter** 
In landing performance measured In terms of touchdown posi- 
tion, sink rate, and speed dispersion. 10 to 15 seconds 
was generally found to be marginal resulting In some scatter 
particularly In turbulent air. 15 to 20 seconds was found 
to be quite satisfactory for all vehicles. Times greater 
than 20 seconds were seldom experienced and cannot be 
evaluated. 

Since the minimum touchdown speed Is determined Inde- 
pendent of the trajectory. Integration of the drag equation 
over the final glide allows the prediction of the speed 
required at flare completion for any desired deceleration 
time. 

Equation (18) In the form 



expresses the final glide deceleration for parabolic drag 
assumptions. 

Equation (19) Integrated Is 



(19) 




where 



<y __ 




( 20 ) 






/ 

m — 
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where 






/^- A-^u. ■^(■^■n) -f ■^■/T u. 









a-f-fzu. 
-^’ V /Z~u^ 



O. =» ~/(7*f 

^ a-Yzou 

Computation with equation (20) Is unhandy and a simple 
graphical solution might be constructed for a particular 
value of E by plotting backwards from touchdown using 
a nominal ^ of about minus 2 degrees. The function Is 

7 

quite linear In the range of Interest and sufficient accu- 
racy can be obtained by plotting only a few points. FIG. 9 
is the approximate solution of equation (20) for the ve- 
hicles tested assuming touchdown at .88 V^. Results of 
the experiment show very close agreement with this analysis 
tending to show the validity of the parabolic drag polar 
assumption, FIG, 10 shows a comparison of trimmed aero- 
dynamic characteristics for vehicles B and C with their 
parabolic approximations. For these examples, the region 
of interest for flare and landing lies between Cj^ of 
approximately .15 and ,40. This region shows quite rea- 
sonable agreement. 
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The deceleration characteristics of re-entry vehicles 
and their landing speeds define the problem In the final 
glide where a rather fixed time must be allowed. Further- 
more, the deceleration Is approximately constant. Consider 
then the form of the drag equation for final glide 

^ ^ ^ ^ -f- 

c/i C U* 




Substituting the velocity ratio at 
(22) / 

I 

Neglecting the gravitational force for simplicity 
( 23 ) ^ 

72 - 



Thus It can be seen that the vehicles with higher lift to 
drag ratios will have superior deceleration characteristics. 
Next consider (23) Integrated between flare completion and 
touchdown. 



(24) 






But If equals .88 Vg 

(25) 







19 



This simplified equation shows the effect of vehicle para- 
meters on final glide and landing speed. The final glide 
characteristics are necessarily determined by both decel- 
eration, a function of E alone, and the speed for minimum 
drag, a function of both wing loading and the drag char- 
acteristics. An example of this is shown in the table 
below. Values shown vjere calculated for T of 15 seconds 
using FIG. 9 » 



VEHICLE 


E 


W 

S 


K 


Cd 

0 


Vr 


Vrain 


Vfg 


A 


3.0 


22.7 


.463 


.06 


230 


202 


317 


B 


5.0 


29.4 


.318 


.0316 


280 


246 


286 


C 


4.4 


25.0 


.381 


.034 


265 


233 


288 


D 


4.4 


31.4 


.258 


.050 


265 


233 


288 


Note: K 


and C 


Do 


based on 


, the 


idealized drag polar. 


The 


above 


table 


depicts 


four of 


the five 


vehicles 


tested 



Vehicle B-, is identical to B except in the computer program- 
ming. Vehicles B and C are representative vehicles of the 

are 

re-entry class vjhile A and ly hypothetical vehicles. 

It may be noted that vehicles C and D, although they 
have different vehicle parameters, show identical landing 
characteristics. This is possible since Vg is defined by 
3 vehicle parameters while E is defined by 2 . Comparing 
vehicles B and C shows the effect of wing loading. Vehicle 
B is superior in maximum lift to drag ratio but must land 
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at the higher speed. Comparing B to A shows the same 
tendency but even more graphically. However, vehicle A 
must be provided the greater starting speed and distance 
from landing area. From the simulator study It could not 
be determined what the effect of the higher deceleration 
and greater landing speeds would be on pilot opinion of 
the vehicle. The measure of success seemed to depend 
primarily on the T, If the vehicle could complete the 
flare with sufficient kinetic energy for approximately 
15 seconds or more In final glide, the landing could be 
successfully completed Independent of the relative decel- 
eration and speed characteristics, 

6 . Flare . 

The flare maneuver (PHASE II) was found to be best 
characterized by the normal load factor used during the 
maneuver. Although It was found that the flare was essen- 
tially a "closed loop" maneuver In which the pilot senses 
his airplane's flight path curvature and altitude, the 
maneuver tends to be accomplished at nearly constant load 
factor. In this experiment, the pilot was not provided 
an accelerometer, yet repeatedly the flare found best liked 
In terms of controllability and accuracy was about 1.2 to 
1,3 "g'*» A series of typical time histories of ^ 
calibrated In "g" units are shown In FIG, 11. 
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Representing "favorable” or near "optimum" pilot 
comments, these data show the consistent average of the 
1.2 to 1.3 "g" flare for 4 out of the 5 vehicles over a 
variety of Initial velocity and altitude conditions. 
Moreover, the simulator tests showed that the favorable 
or most preferred flare load factor was actually Independ- 
ent of altitude and velocity of Initiation, the pilot tend- 
ing to adjust the flare Initiation altitude to give the 
preferred load factor over a number of trial runs, and 
without reference to an accelerometer or the data of that 
particular series of test runs. Flares using higher "g" 
tended to be difficult to control resulting most often 
In misjudged level-off altitudes and subsequent poor land- 
ing or even catastrophic over shooting and premature ground 
contact. No time history data Is available on vehicle A, 
but the same tendency was evident In simulated flight 
experience. 

Analysis of the flare maneuver began by considering 
equation (18) 



( 26 ) 




together with the lift equation (1) In the form 



( 27 ) 
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For 












equation (2?) is 

(28) (07-/J 

u ^ f 



Then 

(29) 



c/<)' r c>7- 0 u. 

C/U^ ( Z£rxU^^U.^-h}f'^ 




To show the effect of the constant load factor assump- 
tion, a series of trajectories were calculated for comparison 
with actual test trajectories in V versus ^ space as re- 
corded on the X-Y plotter (FIG. 12). It may be noted that 
the n equal to 1.2 step function approximation to the actual 
’’closed loop” load factor profile gives excellent results 
In predicting the actual trajectory. 

It would appear that the higher ”g” trajectory would 
be superior since the conservation of kinetic energy for 
use In the final glide is paramount. However, the problem 
of precisely controlling the flight path curvature at the 
higher turning rates, lends favor to the less severe flare. 

Equation (29) reveals that the slope of the trajectory 
in V versus space, which Is related to the very Important 
kinetic energy loss In the flare. Is a strong function of E, 
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To Illustrate this, a family of trajectories for vehicles 
A, B, and C were constructed using the load factor 1.2 ^ 
consistent with the empirical data (FIG. 13). All vehicles 
are started with the steady glide at 450 fps. The final 
glide condition In each case Is minus 2 degrees. Equation 
(29) was used and the non-dlmenslonal u converted to the 
more familiar co-ordinate V. 

It may be observed that the kinetic energy lost during 
the flare Increases rapidly as E decreases. Since the 
steady glide speed Is limited by compressibility to about 
^75 and the flare completion speed is limited by the 

necessity of a finite deceleration time, there Is a defi- 
nite limit to the capability of a vehicle to perform the 
landing maneuver using this technique. This was found In 
test experience with vehicle A wherein the T final glide 
was found to approach zero and landing performance was 
unsatisfactory. Once In the final glide the lower E char- 
acteristics (which might be typical of a vehicle with gear 
and flaps extended) does not seem to be a serious limitation 
If sufficient Initial velocity can be provided. Use of a 
higher load factor In the flare decreases the kinetic energy 
loss, but again the precision required seems to limit this 
technique. 

Another serious vehicle limitation might be the alti- 
tude lost during the flare. If the flare must be Initiated 



24 



at a relatively high altitude, the accuracy with which the 
pilot can perform the maneuver must suffer. In this re- 
spect the simulator, even with Its "Ideal” altimeter and 
limited physical cues, showed this to be the case. As the 
flare altitude Increased, the pull-up maneuver tended to 
become "open-loop" until much of the altitude was dissipated. 
At some lower altitude a second or modified flare was 
usually evident. In general, the high altitude maneuver 
was more difficult and less precise. 

To calculate the altitude loss In a flare consider 
^ 30 ) J/ 

or 

/ \ c/^ _ ^ </^ 

c/t 

with equation ( 28 ) 

( 32 ) ^ . 

This calculation may be carried out stepwise together with 
the trajectory calculation. However, by assuming 
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Then 



(33) 




With equation (33) the altitude lost In the flare for the 
three trajectories of FIG. 12 is 



These results are consistent with the experiment and accurate 
enough for planning an approach since the pilot will actually 
vary his pull-up considerably depending on his Judgment. 

Simulator experience with the flare maneuver through 
more than about 2,000 ft. showed them to be difficult to 
perform precisely. Flare maneuvers at less than about 600 ft. 
also proved to be undesirable but these are seldom encountered 
If the 1.2 "g” pull-up Is used. 

7. Final Approach. 

The steady glide flight conditions of the final approach 
phase (PHASE I) can be determined by calculating a desired 
trajectory backwards from touchdown until an equilibrium 
glide flight condition Is reached for any given vehicle. 

The speed, flight path angle, and rate of descent In 
the final approach will, of course, depend on the lift-drag 



VEHICLE C 



VEHICLE A 



VEHICLE B 



4?80 ft 
1398 ft 
1620 ft 
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characteristics of the vehicle. However, the speed is 
bounded by compressibility effects and the maximum accept- 
able rate of descent probably depends on pilot opinion. 

The latter limitation was not a difficulty in the simulator 
and therefore could not be evaluated. 



8. Synthesis of an Approach. 

An example of how an approach could be synthesized 
is given here and shown in FIG. 14. Vehicle C is used for 
this example. 

It was assumed that the touchdown speed was not other- 
wise limited and .88 Vg was used. It was further assumed 
that 15 seconds was desired in final glide, and the lift- 
drag characteristics were constant (no gear and flaps 
lowered in PHASE III). Final glide angle was minus 2 
degrees. A 1.2 ”g” flare was planned. 

For these conditions 






IS 



.£or 



From FIG. 8 ^ \ 

^ SS ) 



With equation (29) a stepwise trajectory can then be con- 
structed (FIG. 13). The steady glide conditions are found 
to be as follows: 



2 ? 



V.'ith equation (33) 



K- 

K 

_L 

A- 



-fps 

-/r.TS'' 

j./z 



__ // 



This approach was flox-jn on the simulator and the 
actual envelope of 4 trajectories is also shown in FIG. 14. 
The altitude loss ivas 900 ft. and T varied betvjeen approxi- 
mately 13 and 1? seconds. 

With this very simple analysis the ability of any 
particular vehicle to perform a planned approach can be 
quickly evaluated. Simulator experience has shown that 
an actual approach can be very closely approximated by 
this simple method of synthesis. 

The long period dynamics of the re-entry glider land- 
ing can be characterized by eight variables. These are: 

L - Lift V - Velocity 

D - Drag S - Reference area 

g - Acceleration of gravity ^ - Air density 

iZ - Flight path angle W - Vehicle weight 

This study has shown that for typical vehicles the lift- 
drag characteristics can be expressed in terms of a para- 
bolic equation for the landing flight regime. If this 
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Is done» the landing can be characterized by 5 non-dlraenslonal 
groups : 

l. 

c= — =• • — 

JI- 

The flight path angle ^ varies from some steady glide 
value given by 

04 ) / 

5/>r dj- = 

to a typically small final glide angle of about minus 2 
degrees. 

The normal load factor ”n" Is pilot selected and 
essentially constant during the flare maneuver. 

The velocity Is normalized by which contains pure 
vehicle parameters. 

The maximum lift to drag ratio E is a pure vehicle 
characteristic, but Is Independent of wing loading allow- 
ing a degree of freedom between the 2 vehicle parameters 
E and Vjj. 

The deceleration parameter 'T' Is a non-dlmenslonal 
time In which the effects of vehicle parameters on a de- 
sired deceleration time In final glide are shown. 
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FIG. 9 u versus T for final glide, and FIG. I 5 for 
the I Inal approach and flare show how the 5 non-dimensional 
parameters define generalized trajectory characteristics. 
From these generalized characteristics the effect of the 
vehicle parameters E and can be shown as in FIGS. 16 
and 17 . Increasing E gives lovjer final approach glide 
angles and sink rates, less kinetic energy and altitude 
loss in flare, and longer deceleration times in final 
glide for a given initial velocity. Increasing results 
in increased final approach sink rates, greater loss of 
altitude and kinetic energy in the flare, higher landing 
speeds but a slightly Increased deceleration time in final 
glide for a given Initial velocity. 

Since the vehicle parameters E and V are enough to 
separate the trajectory characteristics of vehicles which 
reasonably fit the assumption of parabolic drag polars, it 
seems to suggest that E and could be used as a basis 

II 

for rating vehicles in the manner of flying qualities 
parameters. If a pilot did not object to the short period 
characteristics of a re-entry vehicle, the long period 
dynamics would then become the more important and dominate 
his opinion of any vehicle. Simulator tests with Vehicles 
C and D which have equal E and Vg but varied wing loading 
and drag polar indicated this v/as true. In this case no 
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significant difference In vehicle performance was noted. 
However, the experiment was decidedly limited In this area 
and needs verification. It seems quite reasonable that 
Increases in E would always be welcomed, but that an opti- 
mum In Vjj might be found since the effects of Vg are not 
encountered consistently In all 3 phases of the landing 
as they are with E. 

9 . Landing Dispersion. 

The average landing dispersion for the 4 vehicles 
tested was ^300 feet based on 68 landings In which data 
was available. The maximum deviation from the mean touch- 
down position for each vehicle was as follows: 





Over 


Short 


VEHICLE A 


+1000 


-900 


VEHICLE B 
(Including Bj^) 


+1000 


-1500 


VEHICLE C 


+900 


-1000 


VEHICLE D 


+1000 


-900 



Because the touchdoKTi and flare completion speeds 
were predictable for a given final approach condition. It 
was found that the approach glide could be "aimed” Initially 
to give the vehicle a high probability of touchdown within 
a small landing area (il500 feet from a predicted touchdown 
point). In the presence of turbulence the probability that 
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the airplane might be "bumped" Into an early landing In 
the last 2 or 3 seconds of the final glide was found to 
exist. 

Because of the low altitude, and distance from the 
field (approximately 8,000-12,000 feet) In final glide, 
the pilot may not have a clear view of the landing area. 

For this reason. It Is possible that approach zone range 
markers could be an Important landing aid to Improve touch- 
down accuracy. 

10. Short Period Dynamics. 

Although the primary emphasis In this study was on 
long period vehicle motions. It was found In preliminary 
simulated flights that the basic airplane had unsatisfac- 
tory short period dynamics In the landing flight range. 

The basic re-entry vehicle has of Itself little or no 
rotary damping. Therefore, without some pitch rate feed- 
back, the short period characteristic motion did not pro- 
vide the necessary tight attitude and flight path control 
necessary for precise landing. Particularly In turbulent 
air, the unaugmented longitudinal dynamics proved unsatis- 
factory. Moreover, It seems unlikely that a re-entry 
vehicle configuration could perform Its varied mission 
satisfactorily without extensive longitudinal stability 
augmentation. 
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Using a constant speed approximation, the simplified 
short period motion is given by the characteristic equation 

(35) A* + (v‘- Me- L«.rte)-0 

from REF. 6. The dimensional stability derivatives and 
short period characteristics for Vehicles B and C are 
given In TABLE II for speeds of 350 and 230 fps. Due 
to equipment limitations was Invariant with flight 
condition. 

The damped natural frequency, as approximated, showed 
fairly good agreement with test data. It may be noted 
that the short period characteristics for these vehicles 
without the pitch damper would normally be rated In the 
"poor” to "bad” region In the "thumbprint" (REP. 7)» It 
was not particularly desired to overdamp the short period 
motion but the computer limitations dictated fixed values 
of . 

The damped short period characteristics were found 
satisfactory for the landing study. The minor differences 
between the short period characteristics of the various 
vehicles was negligible so far as the landing problem Is 
concerned. The same was found with the difference in 
control power, gust sensitivity, and Future studies 

might consider these variables in more detail. 
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CONCLUSIONS AND RECOMMENDATIONS 

Tne following general conclusions are made for the 
range of parameters and assumptions of this report: 

1. The final approach and landing technique for the re- 
entry may be divided into 3 phases as follows: 

PHASE I - PINAL APPROACH. Steady state glide at con- 
stant airspeed and flight path angle aimed at a fixed 
geographic point short of the runv;ay, 

PHASE II - FLARE. A nearly constant load factor pull- 
up during which the sink rate is retarded vjhile the 
altitude and flight path angle are reduced to near 
level flight at low altitude, 

PHASE III - FINAL GLIDE. Deceleration to touchdown 
in essentially one "g” flight. 

2. For the landing maneuver, the re-entry vehicle aero- 
dynamics may be approximated by a parabolic drag polar. 

3. For the parabolic drag assumptions, a re-entry vehicle 
approach and landing may be synthesized based on 5 geti- 
eralized, non-dimensional parameters; n, 'if ^ u, E, and T, 
The approach is calculated in reverse from the touchdown 

to the steady state final glide matching boundary conditions 
between phases. 
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4. If the minimum touchdown speed is not otherwise limited, 
a speed of approximately .88 is the minimum landing 
speed giving acceptable final glide performance. 

5. A minimum deceleration time in final glide of 10 to 15 
seconds is required for marginal landing performance while 
15 to 20 seconds yields satisfactory performance. 

6 . An average load factor of 1.2 to I .3 "g" ;vas pilot 
selected as "near optimum” for the flare maneuver independ- 
ent of altitude loss. For trajectory analysis, a step 
function of load factor ’with amplitude of 1.2 "g" gives 
excellent results, 

7 . An average landing dispersion of ^300 feet and a maxi- 
mum of -I 500 feet from the planned touchdown position 
might be expected based on 68 landings with 4 vehicles 
initially positioned on an optimum final glide, 

8 . A maximum sink rate at touchdown of 10 fps is a prac- 
tical limitation and v;as seldom exceeded for an optimum 
approach. Touchdown usually will occur at a small but 
finite flight path angle of about minus 2 degrees. 

9 . The short period dynamics of the basic vehicle are 

generally unsatisfactory for precise landing performance. 

V/ith pitch rate damping incorporated, the long period 

dynamics predominates the landing performance evaluation 

with E and V as the most significant vehicle parameters. 

R 



35 



The follov/ing general recommendations are made as a 
result of this investigation: 

1« Further simulator study to determine the extent to which 
E and determine a re-entry vehicle’s landing performance 

as seen by the pilot; i.e., make v;ide variations of the wing 
loading and parabolic drag equation parameters while keeping 
S and constant. 

2. Further investigation to determine if T final glide is 
a function of E. and Vj^; i.e., study the effect of decelera- 
tion, landing speed, and distance from touchdovjn at flare 
completion. 

3 . Determine if an average flare load factor of 1.2 to 
1.3 "g" is valid over a wide range of pilots and vehicles. 

4. Determine if limits of E and V are defined by approach 

n 

glide sink rate, glide angle and altitude lost in flare. 

5« Investigate the Importance of slope of the lift curve 
on airplane touchdovm attitude and pilot visibility. 
Investigate the short period dynam.ics to find the limit- 
ing pitch damping necessary for precise longitudinal con- 
trol. 

6. Determine the effect of large airspeed and altitude 
offsets from the optimum final glide conditions. 
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•TABLE I 

PROGRAMED VEHICLE CHARACTERISTICS 



QUANTITY 


UNITS 


A 


B 






^ ! 

t 




ft. 2 


160 


136 


136 


160 


160 




lbs. 


3640 


4000 


4000 


4000 1 


5025 

j 




lbs/ ft.^ 


22,7 


29.4 


29.4 


25.0 


0 /• 




P 

slug-ftf 


3698 


4520 


4520 


3698 


3698 ' 




ft. 


21.17 


21.05 


21.05 


21.17 


21.17 






0 


+ .10 


+ .10 


-.12 


-.12 


c. 




2.15 


1.725 


1.72 


2.15 


2.15 






.814 


.573 


.573 


.814 


.814 






.0588 


.0377 


.0564 


.0575 


.0402 






0 


.094 


0 


-.390 


0 






,0227 


-.0306 


.228 


-.158 


-.0349 






1.195 


1.752 


1.21 


1.49 


.391 






0 


2,02 


0 


+ .40 


0 






0 


-.0072 


-.0040 


0 


0 






-.187 


-.0552 


-,0612 


-.187 


1 -.187 






-.268 


-.1245 


-.1146 


- . 268 


-.268 


\ 




0 


-,0524 


0 


0 


0 


/O jc 




- 860 


- 1260 


- 1260 


- 860 


- 860 






V 


V 


V 


V 


V 



NOTE: C ~' = eaulvalent rotary damping for - -10 

UiQ 



TABLE II 



DIMENSIONAL STABILITY DERIVATIVES 



AND 

APPROXIMATE SHORT PERIOD CHARACTERISTICS 





VEHICLE B 


VEHICLE C 






V fps 




V fps 


Derivative 


230 


350 


230 


350 


i 

1 

j 


.527 

-10 

0 

-2.2 


*795 

-10 

0 

-5*1 


.757 

-10 

0 

-10.8 


1 

1.15 ' 

-10 
0 

-25 


I Frequency 
and 

Damping 

Ratio 












.178 


*177 


.115 


.115 




1.^82 


2.25 


3*29 


5.0 


(Ke-iq) 


2.21 


1.50 


1.25 


*923 




2.38 


3.6 


4.3 


6.05 
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APPENDIX A 

Details of the Simulation Equipment 
A,1 Mechanization of Equations. 

The longitudinal equations of motion were simulated 
using the Electronics Associates Pace analog computer and 
a Llbratrol 500 digital computer linked through the analog 
digital voltage converter. The three longitudinal equations 
for body axes plus kinematic relations are: 




(A-3) 




(A-4) 1 / 



(A-5) X y - Y Yo AY 



(A-6) 



(A-7) c/=- ~ AY 



A2 



where the follov;lng approximations have been made: 



2 

^ <h” 2 <Jr" 

Cc::^2> •= ^ A<f 

A2dA<f<2^ 



The sohematic of the mechanization of the above equa- 
tions is shown in FIG. Al. 

The aerodynamic coefficients given in TABLE I vjere 
determined by a best fit to the aerodynamic characteristic 
in the range of both trimmed and untrimmed flight conditio 
suitable to the problem. 

Dynamic pressure was computed in the Libratroi 500 
using analog inputs of velocity and altitude and atmos- 
pheric density programmed in the computer. A standard 
exponential atmosphere of the form 



was used in the computation. 

Experimental data was recorded on the analog's X-Y 
plotter and the 6 channel Visicorder 906C strip chart 
recorder having a galvanometer natural frequency in excess 
of 120 cycles per second. The X-Y plotter was normally 
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employed to plot trajectory data such as: 

^ KS ^ 

The 906 c recorded time histories of the follov;lng variables: 
< 3 ^ A j, ^ a.'^rc/ 

A. 2 Cockpit Display. 

The overall viev: of the cockpit including' fixed base 
seat, instrument display, CRT runway display, and side 
controller is shown in FIG. 1. The detailed sketch of the 
instrument panel and runvjay display is shown in FIG. A2. 

The runvjay display, computed in the Llbratrol 5 OO 
using analog inputs, shov;s a runway 10,000 feet by 200 feet 
in standard perspective. Airplane pitch attitude is also 
an input to the CRT, 

The ^ indicator shows airplane flight path angle in 
degrees. The airspeed indicator was calibrated in feet per 
second. The attitude indicator read in degrees. The 
micter was calibrated to indicate glide path error in hundreds 
of feet. The direction of the error v;as such that the air- 
plane must fly toward the needle. The altim.eter was cali- 
brated to read in thousands of feet to an altitude of 2,000 
feet, then rescaled to hundreds of feet. The roll attitude 
indicator v;as not used in this study. All Indicators were 
scaled for Inputs from the analog computer. 



A, 3 Computation of Runway Display 



The runway display v/as computed in the Libratrcl 5 OO 
using Inputs of h and x from the analog computer. The 
equations of the perspective can be developed considering 
FIG. A3. 

Computing the front and rear of the runway: 










Then 




(A-8) 



(A-9) 




Computing the slope of the runway edges: 

J/ 

~/i2yyt A = =~ 





Then 

(A-10) 
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Using programmed runway dimensions and eye to screen dis- 
tance, the Llbratrol ^00 computed the perspective view 
with equations A-8, A-9, and A-10. An Input of airplane 
attitude ^ to the CRT completed the display. 

A, 4 Simulation of Atmospheric Turbulence. 

The atmospheric turbulence was basically a uniform 
probability distribution of angle of attack with adjustable 
RMS amplitude and bias. The random number capability of 
the Llbratrol ^00 digital computer was used to generate a 
square wave signal F(t) with finite pulses, equal spacing, 
random amplitude and a frequency of 12 cycles per second 
as shown In FIG. A4(a). The probability distribution Is 
shown In FIG. A4(b). From REF. 8 the power spectral den- 



as shown In FIG. A4(c). 

As shown by schematic FIG. A5(a), the random signal 
F(t) was filtered and used to modulate a IKC carrier, then 
recorded on a Sony 262 C AM tape recorder. The corner 
frequency of the first order filter was 1 radian per second. 

The random signal from the tape recorder was then 
demodulated , biased for equal positive and negative and 



slty /^(f) is 
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summed vjith to give the desired turbulence level. The 
gust level was normally adjusted to an equivalent RMS 
vertical gust of 6.5 fps at the speed for minimum drag. 
The RC constant of the demodulator was .5 seconds. The 
filtered and biased turbulence power spectral density and 
probability density are shown in FIG. A6. 
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